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Abstract. Equations of motion are established for a dynamical system in which a 
spacecraft flies close to and interacts with an outer planet and one or more of its 
satellites. For the computation of the state and mass partials (dr.^dr^ and 
dr.^/dni.) needed in the orbit corrections, a set of variational equations is derived. 
These assume a notably compact form through the introduction of a matrix operator 
D. The above system of differential equations is integrated numerically on a 
computer. 


Spacecraft- satellite direction measurements accurate to +10" were simulated 
along three representative Mariner/jupiter/Satum trajectories approaching Io, 

Titan, and Iapetus to within 41, 000, 13,000, and 7, 000 km, respectively. Forex- 
ample, from measurements distributed evenly at half-day intervals over a 60-day arc 
centered on encounter, but none so close that the satellite would fill more than 0. 5° 
in the sky, the orbit of the satellite and that of the spacecraft can be estimated to 
about 100 km. In addition, the mass of the satellite is obtainable to 2.6% for Io, 1.4% 
for Titan, and 9% for Iapetus. If only measurements up to 3 days before satellite en- 
counter are included, the orbit of the satellite or that of the spacecraft can be 
estimated to about 300 km, all information on mass being lost. 
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1. Introduction 


The Mariner/Jupiter/Saturn (MJS) 1977 missions will encounter several natural 
satellites. The goal of this study has been to determine the potential value of space- 
craft-centered TV-imaging data on these satellites for improving the ephemerides of 
the satellites and for aiding the navigation of the spacecraft in their vicinity. More 
specifically, we seek answers to the following questions: 

1) On the basis of simulated satellite imaging data, how large are the improve- 
ments we can expect in the values of (a) the orbital parameters of the satellites and 
of the spacecraft, (b) the masses of the satellites, and (c) the masses and zonal- 
harmonic coefficients of the primaries? 

2) In order to maximize the return of information for question 1), how should 
spacecraft trajectories be chosen and the observations be distributed along them? 

The power of TV-imaging data of satellites to improve ephemerides has already 
been demonstrated (Born and Duxbury, 1972) for Mars' satellites Phobos and Deimos, 
which were photographed from the Mariner 9 spacecraft/orbiter last year. It was 
estimated that whenever, on a TV frame, stars were visible together with one of the 
satellites, the latter's position could be measured to an accuracy of about 5". This 
made it possible to estimate the positions of Mars' satellites to better than 10 km. 
These results are very encouraging, but a direct extrapolation to the satellites of the 
outer planets must be cautioned against. Even if we exclude satellite imaging data 
from the Mariner 9 orbiter phase, which will have no counterpart on a MJS mission, 
there remain significant differences between the two with regard to the satellite data 
and their use: The latter mission is characterized by approach velocities approxi- 
mately 10 times greater; the data must be taken at much greater distances; and the 
target satellites are hundreds of times larger. While a relatively simple analytic 
satellite theory (Aksnes, 1972) was sufficient for analyzing the data on Phobos and 
Deimos, the satellites of Jupiter and Saturn present a much more complicated dynam- 
ical problem owing to mass interactions between them and the spacecraft. This 


2 



circumstance should be welcomed because the deflection of the spacecraft caused by a 
nearby satellite promises to yield a strong determination of its mass. It does mean, 
however, that a numerical integration of the equations of motion and their variational 
equations is the only technique available. Although demanding in computer time, this 
numerical approach is feasible because time intervals of only a few months at each 
planet are involved. 


2. The Dynamical Model 


2. 1 Equations of motion 

For our purpose, we assume that a massless spacecraft (or an orbiter) is moving 
under the attractions of an oblate planet and n-1 gravitationally interacting satellites, 
each treated as a point mass. In Figure 1, m^ denotes the mass of the planet, and m. 
and r. are the mass and the plan eto centric position vector, respectively, of the i-th 
body, i being 1 for the spacecraft and 2, 3, . . . , n for the satellites. The equations of 
motion then take the form 


d 2 ?. 


r. 


dt 


r + (un + n,)4 




1 = 1 , 2 ,.. 


(i) 


where p. denotes the product of the gravitational constant and m., and where r^ . is the 
distance between the i-th and j-th bodies. Furthermore, VjRj denotes the gradient, 
with respect to the cartesian coordinates x^, y., z., of the i-th body, of the disturbing 
function 



P k (sin 6.) . 


(2) 


Here, R will be recognized as the equatorial radius of the planet, as the latitude of 
the i-th body above the planet's equatorial plane, and as the coefficient of the k-th 
zonal harmonic. 
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Let us now assume that the coordinates are referred to the Earth's mean equator 
and equinox of 1950. 0. To express sin 0. in terms of x., y., z., we introduce the right 
ascension ot and the declination 6 of the planet's north pole. We then define a planet 
equatorial system by rotating the xyz-system ot + 90° about the z-axis, and 90° - 6 
about the x-axis. Hence, 



-sin ot , 
-cos a sin 6 , 
cos ot cos 6 , 


cos ot , 
-sin ot sin 6 , 
sin ot cos 6 , 



(3) 


from which transformation it follows that 


r. sin = z/ = cos ot cos 6 


x^ + sin a cos 6 • y^ + sin 6 • z. . 


Then, 


*i R i 


x. /r. , dsin e./dx. 

y. /r. , dsin e./dy. 
z^r. , dsin 


dR i/ dr i 

dR.y'd sin 0. 


(4) 


(5) 


In the Explanatory Supplement to the Ephemeris, the quantities ot and 6 for the planets 
are given as slowly varying linear functions of time. For simplicity, we assume in 
our simulation study that a and 6 are constant during the short planetary-flyby periods, 
and the coordinates will be referred to the equatorial planes of the planets by making 
ot = 0° and 6 = 90°. Equation (5) can then be written 


/2 1 

x. , -x.z./r. 
i ’ li/i 


OO 

, 0 y> + x) Rkp k (sin e i ) / rk+3 

/ 2 


k=2 

y i ’ - y i z i/ r i 


OO 

z. , -z?/4 \ + 1 
1 ’ 1/1 


- ,1 o2 J '' Rkp k (Si " 9 l ) / r i I+2 

1- —* 

k=2 


( 6 ) 
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2. 2 Variational equations 


Let a solution of the equations of motion (1) be expressed symbolically as 


r. 

i 


r ,-( 5 ^*05 ■ ■ • ) Pi J Po) • * • f P m > > 


i = 1,2, ...,n , 


( 7 ) 


where P^. ( j = 1, 2, . . . , n) is the column vector | x^ Q , y^, z. Q , x. Q , y. Q , | defining the 
initial position and velocity of the j-th body, while p^, p 2 , . . . , p m denote m scalar 
parameters, e.g. , satellite masses and zonal-harmonic coefficients. By varying the 
arguments of r., we obtain 


1 

o 

41 

1 


6? 2 

= 

6r 


n 



G ll’ '• *’ G ln 


G nl» * ’ * » G nn 


6P, 


6P, 


6P„ 


n 


g ll> • • g lm 


6p 1 

— — 


6p 2 

g nl’ ’ " * ’ g nm 


1 

Cn 

J* 

i 


or, in matrix form, 


6r — G • 6P + g • 6p 


(8) 


Here, G and g are n x n and n x m matrices whose elements G^ and g^ are, in turn, 
the 3x6 and 3x1 matrices given by 




" dx. 

bx. 

2 

bx. 

i 


■Jx. - 

sv 





ay i 

iy i 

dy i 




55 jo’ 

•" ,di )0 

’ g ik bPy. 


dz. 

J 

bz . 
2 

bz. 

2 


dz. 


5 y jo’ 

••• 9 

30 _ 


_ 5 V 


j 


(9) 


where i and j = 1, 2, . . . , n and k = 1, 2, . . . , m. If we interpret 5r as the difference 
between observed and computed values of r, then (8) becomes an equation of condition. 
From a series of observations, corrections 6P and 6p to the nominal values of the 
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constant parameters P and p can be found by the method of least squares, provided 
that the coefficients G and g can be computed. To this end, we differentiate Equations 
(9) twice with respect to t: 


• • n 

hr ^-^dr dr 
3 1=1 z 3 


dr^dr. 

g ik ’ 


where the terms under the summation signs are due to the implicit dependence of r. in 
(1) on P. and p^, while the last term arises because of the explicit dependence on p^. 
For the determination of G- and g.^, we thus arrive at the differential equations 
(variational equations), 


n x— 

r — A dr. 


Z or i 

5fT G *j ’ 


j I? 2, • . . y n ? 


.. y' d7 i a7 i 
g ik = Z^S?7 g *k + 55r 


k 1, 2, . . . , n , 


with initial conditions, 


1 0 0 0 0 0 

G„= 0 1 0 0 0 0 6 G.. 

0 0 1 0 0 0 


0 0 0 1 0 0 

0 0 0 0 1 0 6 .. ; 

i] 

0 0 0 0 0 1 


«ik = g ik = 0 > 

0 


1 for i = j 
0 for i ^ j 
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To derive expressions for the partial derivatives in Equations (10), we vary both 
sides of Equation (1): 


6r. =-6(p ( 


i 3 \ )i 3/ 


+ (^0 + l*j) D(r.) 6r. 


D(r j ) 6r i D(r j “ r i> 6(r j “ r i ) +V i 6R i 


( 12 ) 


J 


3 


where D is a matrix operator defined in terms of a vector r = [ x, y, z] and the identity 
matrix I as follows: 


D(r)=-^ 

r 


3rr 


- I 


x 2 - r 2 /3 , xy , 


xz 


2 2 /o 

yx , y - r /3 , yz 

2 2 /o 

zx , zy , z - r /3 


(13) 


In Equation (12), is to be varied only with respect to the quantity in Equation 

(2). Thus, we neglect terms of the type J^Sr., which should contribute only negli- 
gibly during the short time intervals considered here. The terms factored by the 
central mass m-q will of course dominate, except that if r.. should become very small 
during a close approach to a satellite, the corresponding term will be brought into 

sudden prominence. We now immediately deduce 

• • 

d7 i - - - 

5 =-^jD(r.)-^D(r.-r i ) , 


• • 
* r i 

a? 

o 

3. 

II 

+ n i ) D(r.) 

n 

+1} 

i 



i=i 





• • 

dr i 

• • 

ar. 

r. 

j 

• « 

<*i 

a^ 0 

II 

it 

= ~~3 ’ 

r i 



(14) 


3 3 ’ 

r.. r. 

3i 3 


(continued) 
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'Vi/ r ? 

- y i z i/ r f 

-m ♦ 


(k + 1) P k (sin e.)/r^ +3 
-P '(Sin e.)/ r k +2 


where i, j ( i ^ j ) = 1, 2, . . . , n, k = 2, 3, . . . , <= ; P k , the Legendre polynomial of order 
k, should not be confused with the parameters P^ introduced in Equation (7). 


3 ^ Nume rical Resul ts^ 

Based on the equations established in the preceding section, a computer program 
has been developed that can differentially correct the following 17 parameters: 

P 1 = t x 10’ y 10’ z 10’ x 10’ y 10’ z 10^ ’ initial conditions of spacecraft; 

P 2 = [ x 20’ y 20’ z 20’ x 20> y 20’ z 20^ » initial conditions of satellite; 
p^ = Pg , mass of satellite; 

Pg = M-p , mass of primary; 

P 3’ p 4’ p 5 = 3 2 ’ 3 3* J 4 ’ form factors of Primary. 

For simplicity, we have assumed that the satellite under consideration is body num- 
ber 2. Although we have limited the number of bodies to three (primary, satellite, 
and spacecraft) in our simulation runs, the program will accept an arbitrary number 
of perturbing bodies, which need not be only satellites. If we were dealing with 
actual imaging data, such additional bodies would, of course, have to be included. 
However, they would only complicate our simulation without adding much new insight. 

Our equation of condition (8) now becomes 
2 _5 

67 =7} G 2i ■ G li> 8P i + IZ (g 2k - *lk> 6 Pk > < 15 > 

i=l k=l 
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where the G's, g's, and 8r, i.e., the difference between the "observed" and the com- 
puted values of r = r 2 - r^, are obtained by numerical integration of Equations (1) and 
(10) by means of (14). Note that since the spacecraft has no effect on the satellite, 

d7 2 

say, longitude cp and latitude 0 in the planet-equatorial xyz-system. Equation (15) 
must then be modified to 


/ 


dr^ = = 0. In practice, we will be dealing with direction measurements — 


cos 96cp 


-sin cp , 

cos cp , 

0 

6r 

60 


-sin 0 cos cp , 

-sin 9 sin cp , 

cos 0 

r 


(16) 


We have used this program to analyze three representative MJS trajectories 
(Burke et al. . 1972) approaching Io, Titan, and Iapetus to within about 41, 000 
13,000, and 7,000 km, respectively. The corresponding initial conditions, in the 
form of nominal values of the above-mentioned 17 parameters, were used to generate 
fictitious spacecraft/satellite direction measurements with the distributions shown in 
Table I. The table also gives the times of encounter with the planets (E) and with the 
satellites, and the associated minimum ranges (A) to the center of each body. Finally, 
for each case, the range and time of the closest observations are given. The observa- 
tions are distributed evenly within the indicated time intervals, except for a few 
observations dropped because of the imposition of a cutoff value on the satellite range. 
Note that in the first two cases of each flyby, the observations are grouped fairly 
symmetrically about the satellite encounters. As will be demonstrated later, it is 
also important to observe the satellites after encounter. They will then appear much 
fainter than before encounter because of the large phase angles ( > 100°) under which 
they will be seen. Since it is difficult to record very bright satellites and faint refer- 
ence stars simultaneously, the reduced satellite brightness after encounter should 
actually be an advantage. 

It is perhaps unrealistic to hope to obtain precise astrometric data of the satel- 
lites when they fill more than 0. 5° of the field of view of the camera. This occurs at 
distances below roughly 400, 000, 560, 000, and 150, 000 km from Io, Titan, and 
Iapetus, respectively. For this reason, and also because of the high observation 
frequency near encounter (~2 observations/hr), Cases 1 for the former two satellites 
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are probably not realistic. Case 1 for Iapetus may be more indicative of what can be 
achieved in practice. There, observations for which the angular diameter of the satel- 
lite exceeds 0. 5° are excluded, the shortest observation interval being 3 hr. The 
observations are now spread out over 60 instead of 10 days, although one-half of the 
observations are concentrated in 1 week around encounter. In Cases 2 for all three 
satellites, the imaging requirements are relaxed even further by distributing the 
observations evenly over the entire 60-day arc at half-day intervals. 

Cases 3 and 4 should be representative for the imaging data that would be avail- 
able for navigating the spacecraft and aiding camera pointing during the satellite 
flybys. Excluded from Cases 3 and 4 are observations obtained later than about 1 
and 3 days, respectively, before the satellite encounters in order to allow time for 
transmitting the observations to Earth, performing the orbit corrections, and trans- 
mitting the necessary commands back to the spacecraft. 

To simulate noise, a random error within the limits ± 10”, which is compatible 
with the capability of the baseline MJS imaging system, was added to each observa- 
tion. Finally, before the program for each of the cases listed in Table II was run, 
the parameters to be corrected were contaminated with errors of 1000 km and 0.01 
km/s in the position and velocity components, and 10% in the other parameters. 

Table n gives for the corrected parameters the standard errors resulting from the 
least-squares fits, together with their actual errors, i. e. , the deviations from the 
nominal values. It is reassuring to observe that these two types of errors are in 
reasonable agreement. Interestingly enough, probable errors, defined as 0. 67 times 
the standard errors, would fit the actual errors considerably better in most cases. 

To the extent that the observation errors are truely random, the standard errors 
should be proportional to the average observation error (this was actually verified by 
rerunning Case 1 for Io with observation errors between + 20” instead of ± 10”) and 
to the inverse square root of the number of observations. 
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4. Conclusions 


On the basis of Table n and a few cases not reported there, the following answers 
to the questions raised in the Introduction have evolved: 

la) The MJS satellite imaging data ( Cases 1 and 2) should make it possible to 
estimate, for a period of a few months, the positions of Io, Titan, Iapetus, and the 
spacecraft to within several tens to several hundreds of kilometers, depending on the 
distribution of the data and on the satellite. The corresponding velocity errors vary 
from a fraction of a meter per second to a few meters per second. Comparable re- 
sults can be expected for other satellites that the spacecraft approaches to within, 
say, one million kilometers. Even better results are possible by utilizing the radio 
tracking data of the spacecraft in addition to the imaging data. If the observations 
after encounter are excluded (Cases 3 and 4), the position components of the satellite 
become rather strongly correlated with those of the spacecraft. However, if the 
latter parameters are assumed to be known well enough from the radio tracking data, 
the imaging data suffice for improving the ephemeris of the satellite to a few hundred 
kilometers (Cases 3c and 4). This is probably sufficient for navigation purposes. 

The seriousness of the problem of camera pointing is illustrated rather dra- 
matically during the close passage of, for example, Titan in Case 1: At the epoch 
(Saturn encounter), the induced position and velocity errors of the spacecraft relative 
to Titan amount to a few thousand kilometers and a few tens of meters per second. 
These errors are rather smaller than what can be expected in practice. Nevertheless, 
they cause an impressive 17° displacement of Titan as seen from the spacecraft at a 
distance of 14, 000 km. After orbit correction, which in this case needed many itera- 
tions before converging, this displacement was reduced to 0. 003° ! Although, as 
explained earlier, the distribution of the observations is undoubtedly too optimistic in 
this case, it does illustrate to what extent on-board optical measurements can reduce 
the a posteriori uncertainties in camera pointing. 

lb) It may be possible to reduce by as much as a factor of 10 ( Cases 1 and 2) the 
current uncertainties in satellite mass, often quoted as 10% for Io and Titan and 

50 - 75% for Iapetus (de Sitter, 1931; Jeffreys, 1953; Kozai, 1957). Since the 
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satellite mass is estimated from the perturbing effect on the spacecraft, the mass 
estimate will depend rather strongly on the flyby distance to the satellite. 

lc) An attempt to solve for the mass of Jupiter failed because of very strong 
correlations with the mass and velocity of Io. While J 2 of Jupiter (Case lb) could be 
determined to + 14%, this cannot compete with the 1% uncertainty that de Sitter (1931) 
obtained by analyzing the motions of the Galilean satellites. The addition of J g and J^ 
to Case lb gave uncertainties of almost 100% in these parameters. In view of these 
negative results, no attempts were made to solve for the Pq. Jg, Jg, and J 4 of Saturn. 

2) It is highly desirable to extend the observations past the encounters with the 
satellite and its primary. Not only will this increase the total length of observing 
time, during which the satellite orbit will be seen under widely different viewing 
angles; but also it becomes possible to exploit the pronounced bending of the space- 
craft's trajectory during the two encounters. We have seen that the bending caused by 
the satellite can reveal a surprisingly accurate value of its mass. Although a close 
passage of the satellite is required for a strong determination of mass, it is not man- 
datory that the satellite be observed during close-up (compare Cases 1 and 2 for Titan). 
For the improvement of satellite ephemerides, it is not even important that the space- 
craft should come very close to the satellite. As long as their minimum range is not 
decreased (Cases 1 and 2 for Iapetus), it appears that relatively little is gained by 
concentrating the observations around the time of encounter. On the other hand, it is 

probably not worthwhile to include observations that are much farther removed than 

d 

30 days from encounter. At E + 30 , an imaging resolution of 10” would correspond 
to Earth-based resolutions of 0. 5” at Jupiter and 0. 25” at Saturn. 
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Table I 

Distribution of Observations 


Flyby 


Number of 

Observation 

observations 

interval 


Comments 


Jupiter /Io (JST 77) 



Saturn/Titan (JST 77) 



14 








Standard Errors (line 1) and Actual Errors (line 2) 
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Fig. 1. MJS flyby geometry. 
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